A nanomechanical resonator coupled linearly via its momentum 

to a quantum point contact 

Latchezar L. Benatov and Miles P. Blencowe 

Department of Physics and Astronomy, 
Dartmouth College, Hanover, New Hampshire 03755, USA 
(Dated: April 18, 2012) 

Abstract 

We use a Born-Markov approximated master equation approach to study the symmetrized- 
in-frequency current noise spectrum and the oscillator steady state of a nanoelectromechanical 
system where a nanoscale resonator is coupled linearly via its momentum to a quantum point 
contact (QPC). Our current noise spectra exhibit clear signatures of the quantum correlations 
between the QPC current and the back-action force on the oscillator at a value of the relative 
tunneling phase {rj = —n/2) where such correlations are expected to be maximized. We also show 
that the steady state of the oscillator obeys a classical Fokker-Planck equation, but can experience 
thermomechanical noise squeezing in the presence of a momentum-coupled detector bath and a 
position-coupled environmental bath. Besides, the full master equation clearly shows that half 
of the detector back-action is correlated with electron tunneling, indicating a departure from the 
model of the detector as an effective bath and suggesting that a future calculation valid at lower 
bias voltage, stronger tunneling and/or stronger coupling might reveal interesting quantum effects 
in the oscillator dynamics. 

PACS numbers: 85.85.-|-j,72.70.-Fm 



I. INTRODUCTION 



Nanoelectromechanical systems (NEMS), in which a nano-to-micrometer scale mechani- 
cal resonator is coupled to an electronic device of similar dimensions, have received a great 
deal of theoretical and experimental attention in recent years, as these systems are a promis- 
ing tool for gaining a deeper understanding of the quantum-to-classical transition in physics, 
in addition to their useful applications in ultrasensitive metrology.^ A wide variety of NEMS 
have already been realized experimentally, such as, for example, a doubly-clamped nanobeam 
coupled to a superconducting single-electron transistor (SSET),'2a suspended carbon nan- 
otube coupled to an embedded quantum dolPor single-electron transistor (SET),l^a doubly- 
clamped beam coupled to an external SET,'^ and a micromechanical cantilever coupled to 
a quantum point contact (QPC).^ There have also been a number of theoretical studies of 
NEMS, in which the oscillator is coupled linearly via its position to a QPCPISl 

or a SETPit 

has been shown both theoretically and experimentally that the effect of the electronic device 
(detector) on the oscillator is very similar to that of a thermal bath with a certain effective 
temperature and damping constant, even though the detector is in a far-from-equilibrium 
state.'^ Besides, the oscillator can also have a strong effect on the detector, producing a 
Fano-like current noise spectrumP 

Most of the studies conducted so far have focused on a position-dependent linear coupling 
between the oscillator and the detector. It is interesting to see how a momentum-dependent 
coupling changes the oscillator steady state and the detector current noise spectrum, and 
whether there is a nontrivial interplay between the effects of the momentum-coupled ther- 
mal bath associated with the detector and those of the position-coupled bath due to the 
environment of the oscillator. Normally, when a position-coupled detector acts as a thermal 
bath with effective temperature T^et and damping constant 7det5 in addition to the environ- 
mental bath temperature T and damping 70, the oscillator is in a thermal state with effective 
damping 7efr = 70 + 7det and temperature Tes = (70^ + 7det^det)/7efr-'^ In the absence of 
an environment, a momentum-coupled detector is equivalent to a position-coupled one un- 
der the canonical transformation, which interchanges the oscillator position and momentum 
coordinates. However, the unavoidable presence of a position-coupled environmental bath 
breaks this symmetry, leading to potentially new and interesting physics. 

One example of a NEMS which, after an appropriate transformation (see Sec. below), 
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can be described by a momentum-coupled effective Hamiltonian, was studied experimen- 
tally by Stettenheim et al.'^ Their experiment involved a nanomechanical GaAs oscillator 
coupled piezoelectrically to a radio-frequency QPC embedded in it. Measurements of the 
current noise through the QPC detector showed that the quantum statistical fluctuations of 
tunneling electrons could affect the macroscopic dynamics of the host crystal. Fig. [TJ^a) (re- 
produced from Ref. [T6j) shows the GaAs crystal containing a two-dimensional electron gas 
(2DEG). A displacement dy of the front and back faces of the crystal leads to a compression 
dz at the midpoint of the left face and a corresponding expansion dz at the midpoint of the 
right face, as shown. The resulting strain Sy^ = 2dz/w, where w is the width of the crystal, 
produces, through the piezoelectric coupling constant 6x4, a bulk polarization = Cx^Syz, 
which is assumed to be in the direction of transport through the QPC. The 2DEG electrons 
will try to screen the polarization charge, but under the gates and in the QPC, where the 
2DEG is depleted, there will be a net electric field and a corresponding potential difference 
de = Xdz between the left (L) and right {R) reservoirs, leading to a current / through the 
QPC. One of the normal vibrational modes of GaAs has a polarization field as in Fig. [l]^a), 
and thus the QPC current / can provide information on the displacement dz of the crystal as 
it oscillates in this mode. On the other hand, the unavoidable shot noise due to partitioning 
of electron-hole pairs at the QPC leads to charge fluctuations dn in reservoirs L and R, and 
a corresponding back-action force dF = rjdn on the oscillator via the piezoelectric effect, 
completing a feedback loop between the mechanical and electronic degrees of freedom. Thus 
one expects both the mechanical motion of the resonator and the current noise through the 
QPC detector to be peaked at the oscillator frequency. An interesting result in the experi- 
ment, which we set out to investigate theoretically in the present study, is that the current 
noise spectrum of the QPC displays super-Poissonian values close to the oscillator frequency 
yet sub-Poissonian values away from it, indicating bunching and anti-bunching of electron 
tunneling events due to the coupling to the oscillator. One important caveat to keep in mind 
when comparing theory and experiment, however, is that the experiment was performed in 
the strong tunneling regime, where the QPC conductance Gqpc ~ O.SGq and Go = 2e^//i is 
the conductance quantum, whereas our theoretical calculation is based on the assumption 
of weak tunneling. The case of strong tunneling will be considered in a future publication. 

The present paper is structured as follows: In Sec. [IT], we perform a polaron-like trans- 
formation on the Hamiltonian of the system described above, leading to an oscillator 
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FIG. 1. (a) Displacement and associated polarization of the GaAs crystal (Fig. lb in Ref. [16]). 
(b) Schematic diagram of the different electronic reservoirs and tunnel barriers in the system (Fig. 
Ic in Ref. |16j). The symbols are defined in the text. 



momentum-dependent tunneling amplitude across the QPC, and derive a Born-Markov 

we 



approximated master equation for the reduced oscillator density matrix. In Sec. Ill 



solve the master equation to obtain the average current and symmetrized-in-frequency cur- 
rent noise spectrum through the QPC. The noise spectra thus obtained are investigated in 



Sec. IV for a wide range of system parameters. In Sec. M, we use the Wigner representation 



of quantum mechanics to study the steady state of the oscillator. We present our conclu- 



sions in Sec. VI Details of the polaron transformation, the derivation of the Born-Markov 
master equation, and its solution to obtain the current noise spectrum are presented in 
Appendices |A} |B] and [C| respectively. 



II. DERIVATION OF THE BORN-MARKOV MASTER EQUATION FOR THE 
FOURIER-TRANSFORMED REDUCED DENSITY MATRIX OF THE OSCILLA- 
TOR 

As explained in the introduction, the piezoelectric coupling between the 2DEG reservoirs 
and the flexing GaAs crystal can be modeled by a linear dependence of the single electron 
energy levels in the L and R reservoirs on the crystal displacement (Fig. [l](b)). The energy 
levels of the emitter {E) and collector (C) reservoirs, which represent the leads connected 
via tunnel barriers to the L and R reservoirs, respectively, are assumed to be fixed. Thus 
one starts with the Hamiltonian H = Hgyg + i/bath + -f^int, where the system, interaction, 
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and bath Hamiltonians are, respectively: 

Hsys = hwraa^a + Y^{eL - Xx/2)b[bL + J^i^R + Ax/2)6^6r, 

L R 
E C 

Hint = Yl ^^ELblbL + Yl ^^CRb^bRY^ + Yl ^^L^bn + H.c. (1) 

E,L/ C^R LjR 

Here the operators 6j {b\) denote the annihilation (creation) operators for the energy levels 
ei of a given reservoir i = E,L, R, C, the Y operator counts electrons traversing the RC 
reservoir barrier, the Qij are the tunneling amplitudes between the various adjacent reser- 
voirs, the parameter A describes the piezoelectric coupling between the (bosonic) crystal 
vibrational mode x = Xzp{a + a"'') with frequency and the L, R reservoir electrons, and 
is the vibrational amplitude zero-point uncertainty. It is important to note that we use 
X, and not z as in Ref. (TB] and Fig. [l| as the direction of motion of the oscillator throughout 
the rest of this paper. 

The coupling described above is rather unusual - in most electromechanical systems stud- 
ied so far, the oscillator position affects either the charge state of a single island,'^^^^'^ or 
the tunnel-barrier potential,'^^ rather than the reservoirs to the left and right of the tunnel 
barrier. To obtain the current noise spectrum, one could derive a Born-Markov approxi- 
mated master equation directly from the above Hamiltonian, tracing over the bath degrees 
of freedom comprising the E and C reservoir electrons. However, this approach can be quite 
involved due to the large number of coupled second order moment equations that one needs 
to solve. An easier method is to perform a polaron-like transformation on Eq. ([T|, and derive 
a much simpler effective Hamiltonian for our system, in which the coupling maps effectively 
onto a momentwrn-dependent tunnel barrier potential. We replace H — UHW, where the 
unitary operator is 

Ax. 
2hw. 

Expanding to first order in the oscillator displacement and neglecting the quartic terms 
in the bi and bn operators as well as the momentum dependence in the E and C contact 
resistance barrier terms, we arrive at the much simpler Hamiltonian H = ifosc + -f^bath + -f^int, 
where 



U = exp 



(2) 



L R 

Ax. 



1_ _pP(at_a) 

hUm. 



J2^^LRb[bRY + }i.c., (3) 



L,R 

and y now counts electrons tunneling through the LR reservoir barrier. This calculation is 
described in detail in Appendix |Aj In the above derivation, we have assumed weak coupling 
between the oscillator and QPC, A ^ 1, where the dimensionless coupling parameter is 

and m is the oscillator mass, as well as weak tunneling, to ^ 1? where the bare tunneling 
amplitude to is defined below in Eq. ([s]). Together with the assumption of high bias voltage 
across the QPC, eVjhbJm ^ 1, needed to make the Born-Markov approximation described 
below, these are the three main conditions of validity of our calculation. As already men- 
tioned, the experiment of Ref. [TB] is in the regime of strong tunneling, where the polaronic 
and Born-Markov approximations are no longer valid and scattering matrix methods can be 
used insteacP^ - an approach we intend to investigate in future work. 

It is convenient to express ifmt in terms of the oscillator momentum, p = {ih/2xzp){a^ — a)'- 

i^int = T(p)^6jj6z. + H.c., (5) 

L,R 

where 

and we have assumed that VLlr is level- independent. From this point, we follow the approach 
of Doiron,^° since our Hamiltonian has exactly the same form as his, except that his tunneling 
amplitude T{x) is position-dependent. We can write 

f{p) = ^{t, + e^H,p)Y\ (7) 

where A is the constant density of states in the reservoirs, and 

to = 2'KKh\nLRl 

27rA|ilLij|A 



t 
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We have taken the absolute value of fi^^ since the overall phase of to and ti is unimportant 
and only the relative phase difference t] matters physically. In Ref. |H] and most other studies 
so far, only the case of 77 = was investigated for the position-coupled system, implying 
a zero average back-action force on the oscillator. In our case, rj = —7c/2 implies that the 



average back-action force is non-zero (cf. Eq. (13) below; under the canonical transformation, 



which interchanges the oscillator position and momentum coordinates, Fo{r]) indeed becomes 



the average back-action force, cf. the first line of Eq. (10) below). The case of non-zero rj is 
considered in Ref. [TT], where two tunnel junctions, one of which is linearly coupled to an 
oscillator via its position, are arranged in an Aharonov-Bohm-type setup, and the magnetic 
flux through the loop can be used to tune the phase r] between the oscillator-independent 
and oscillator-dependent total tunneling amplitudes. It is shown that when 77 = mod vr, 
the current noise spectrum of the detector is proportional to the position spectrum of the 
oscillator as in Ref. p], but when t] = tt/2 mod vr, the noise spectrum is proportional to the 
momentum spectrum of the oscillator. On the other hand, in Ref. [T^ it is demonstrated 
that for a non- stationary oscillator coupled to a single QPC via its position, the current 
noise spectrum of the detector is complex-valued and contains information about both the 
oscillator position and the oscillator momentum, even when rj = 0. 

Assuming that the oscillator-bath coupling is weak and the bath correlations decay much 
faster than the characteristic timescale of the oscillator, we can use a Born-Markov approxi- 
mation technique to derive a master equation for the Fourier-transformed reduced oscillator 
density matrix 

p(x;t) = $^e^^^p(Ar;t), (9) 

N 

where p{N]t) = {N\posc{t)\N) is the A^-resolved oscillator density matrix and N is the 
number of electrons that have tunneled from the left into the right lead at time t. The 
details of this calculation are presented in Appendix [B} and the final result is 

(T=±l 

+ l^^^li 7a[]3, {x,p(x;t)}] 

a=±l 

- [x,p{x;t)]] - ^7o[£, {p,p(x;^)}] 



(T=±l ^ 



-^(— ^)(iotib,p(x;t)]) 

cr=±l 1 
cr=±l 1 

jr^ [—^) (^o^i{£, P(X; i)}) 



a=±l - ^1 



+ E ^^^^^ ic^-r^) - (10) 



We have defined 



7<T 



( — ) [r^(;ia;m) - r^(-;ia;m)], (11) 
... ^^0^ 

= ^{^j^\T„{nuj^) + T,{-nLo^)l (12) 

Fo(r/) = nsinr7(^)5^ar.(0), (13) 

where cr = ±1 and the forward (left to right) and backward (right to left) tunneling rates 
are, respectively, 

hT+{E)= / de\to\'f{e-fiL)[l-f{e-fiR + E)], (14) 
Jo 

POO 

hT4E)= / de\to\'f{e-fXn)[l-fi^-f^L + E)], (15) 
Jo 

where /ij and fi are the chemical potential and Fermi function of reservoir i. We have also 
modeled the environment of the oscillator as a thermal bath by including external diffusion 



and damping terms (the first and second terms on the third line of Eq. (10), respectively), 
where 

hijj„ 

2m%kBT (when ksT ^ hi^m), (16) 



Do = mjohum coth 



7o is the external oscillator damping constant (defined in such a way that ~2^oP is the 
classical external damping force), and T is the temperature of the environment. 

If the external damping and diffusion terms were not present, our system would be identi- 
cal under the canonical transformation {p -H- mUmx) to the position-coupled resonator-QPC 
systems studied by other groups.'^'^, but with non-zero 77. However, the presence of such 
terms destroys this correspondence and creates a fundamentally new situation, in which there 
is a potentially interesting interplay between the effects of the position-coupled environment 
and those of the momentum-coupled detector on the oscillator. 
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III. SOLVING THE MASTER EQUATION TO FIND THE AVERAGE CURRENT 
AND SYMMETRIZED-IN-FREQUENCY CURRENT NOISE SPECTRUM 



We assume a large forward bias (so that we can drop the a = —1 terms in Eq. (10)) and 
zero temperature in the leads so that we can set /(e — Hl,r) = Q{fJ'L,R — e) in the definitions 
of T±(E) and compute the resulting simple integrals in Eqs. ([I4))-([I5|). The average current 
and current noise can be computed from the moments of N using the formula 

(^'^W)=^-"Tr(|^p(x;t))^^^ (17) 



and taking the x-derivative and trace of Eq. (10) (cf. Eq. (|9])). After a straightforward 
calculation we find for the average current 

(I) = e^(iV) = '-^{tl + tlip^)) - -^^ f^lAx) - em^col^^, (18) 

where 7+ = fi^t\/2mh and eV = fii — f^R is the QPC bias voltage. To find the symmetrized- 
in-frequency current noise spectrum, 

Sriu) = \ l\siiT),5m})e^--dT, (19) 

where 61 = I — {I), we use the MacDonald formula!^ 

Si{u) = 2e^u dtsm{ut)-^{{N\t))). (20) 
Jo 

The time derivative of the variance of can be computed from the time derivatives of the 
moments using the expression 

|((A^'W)) = J^{N'{t)) - 2(iV(t))|(iV(t)). (21) 

One obtains the following result: 

Sj = 2e(J) + A^,, (22) 

where the first term is the Poissonian (oscillator-independent) part of the noise, and the 
oscillator-dependent part is given by the integral 

ASj = 2e^io rdtsm{iot)(- ^^'^^^+ ^((xiV)) + '^tl{{p^ N))) . (23) 



h ti h 

It is possible to solve analytically for the time dependence of the cumulants {{xN)) and 
{{p'^N)), and integrate them to obtain an algebraic expression for ASi. (Note: The double 
angular brackets used throughout this article denote second cumulants, i.e. covariances, 
between products of powers of the oscillator coordinates (e.g. x, p^ or xp) and A^, not higher 
order cumulants.) The full calculation and results are presented in Appendix [Cj 



IV. RESULTS 



To plot the current noise spectrum, it is useful to put the equations in Appendix [C] in 
dimensionless form. We define 

X 

x= — , 

P 

P= , 

Pzp 

u=— (24) 

to be the dimensionless oscillator position, oscillator momentum, time and oscillator fre- 
quency, respectively, where 



Xy, 



2mujr, 



Pzp = y (25) 

are the oscillator position and momentum zero-point uncertainties. Our system is then 
governed by five dimensionless parameters: 



f, _ 7o 
J- — — , 

T = (26) 

and to, i.e. the dimensionless bias voltage, coupling, external damping, external temperature, 
and bare tunneling amplitude. 

We plot the non-Poissonian part of the symmetrized-in-frequency current noise spectrum 
in units of the average current, i.e. ASj{uj)/2e{I) , versus dimensionless frequency u for 
different values of the dimensionless parameters. First, we explore the regime of high bias 
voltage and high external temperature. Each of Figs. [2]-|6]shows plots of the non-Poissonian 
current noise for different values of a certain dimensionless parameter, the other parameters 
being held fixed. The parameters being varied are the bare tunneling amplitude to (Fig. |2]), 
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FIG. 2. Dimensionless non-Poissonian current noise spectrum for Iq = 0.1, 0.2 and 0.3. The 
remaining parameters are V = 2 x 10^, A = 1 x 10^'^, Fq = 5 x 10~^, and T = 1 x 10^. 




FIG. 3. Dimensionless non-Poissonian current noise spectrum for F = 2 x 10^, 2 x 10^ and 2 x 10^. 
The remaining parameters are to = 0.2, A = 1 x 10"'^, Fq = 5 x 10^^, and T = 1 x 10^. 
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FIG. 4. Dimensionless non-Poissonian current noise spectrum for A = 5 x 10 ^, 1 x 10 and 
2 X 10^^. The remaining parameters are to = 0.2, V = 2 x 10^, f o = 5 x 10^^, and T = 1 x 10'^. 




FIG. 5. Dimensionless non-Poissonian current noise spectrum for T = 5 x 10^, 1 x 10"^ and 2 x 10^. 
The remaining parameters are to = 0.2, V = 2 x 10^, A = 1 x 10^'^, and f o = 5 x 10^^. 
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FIG. 6. Dimensionless non-Poissonian current noise spectrum for f o = 0, 5 x 10 ^ and 5 x 10 ^ . 
The remaining parameters are to = 0.2, V = 2 x 10^, A = 1 x 10^'^, and T = 1 x 10^. 



the bias voltage V (Fig. |3j), the couphng A (Fig. |4|), the external temperature T (Fig. |5j), 
and the external oscillator damping f q (Fig- 16|. 

In all the plots, one observes three peaks in the noise spectrum at cD = 0, ±2, as well 
as two resonance-antiresonance features at a; = ±1 (the noise is a symmetric function of u 
and is not plotted for cD < 0). The peaks, especially the ones at a; = 0, ±2, are very sharp 
and can be as high as 10^ for some parameter values, so it was necessary to truncate them 
in order to resolve the off-peak behavior. To give a general sense of the peak magnitudes 
and their variation, the cD = peak ranges from 6 x 10^ to 8 x 10^, the cD = ±2 peaks are 
usually half as high and range from 3 x 10^ to 4 x 10^, the u = ±1 resonance peaks range 
from 2 X 10^^ to 1 x 10^, and the antiresonance features range from —1 x 10~^ to —2 x 10^^. 
The peaks and resonance-antiresonance features tend to broaden and become higher/deeper 
as to, V and A increase (Figs. |2]-[4|. The peaks at a) = 0, ±2 broaden and become higher 
as T increases, but broaden and become lower as Fq increases. At u = ±1, the resonance 
peak gets higher as T increases and lower as Fq increases, whereas the antiresonance gets 
shallower in both cases, and there is no noticeable change in the width of either feature 
(Figs. Island [6]). 
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Our noise spectrum is similar to that derived in Ref. [8] for a position-coupled oscillator 
and QPC with rj = 0, except that we see a resonance-antiresonance feature instead of a 
positive peak at a) = ±1. In the rj = case, the noise spectrum near Co = ±1 is pro- 
portional to the position spectrum of the oscillator.'^ It can be expressed as a leading term, 
corresponding to a classically fluctuating junction conductance, minus a quantum correction 
term, which arises from the correlations between the intrinsic shot noise of the detector and 
the back-action force on the oscillator. This quantum correction is always smaller than the 
leading term when r/ = 0, resulting in a positive Lorentzian peak at a; = ±1. In Ref. [5], the 
non-Gaussian correlations between the junction current and back-action force are derived 
using a simple model, in which tunneling electrons impart random momentum kicks to the 
oscillator at the exact moment of tunneling, the typical size of the kicks being set by the 
Heisenberg uncertainty principle. In the r] = case, the most notable effects of these correla- 
tions are an enhancement of the a) = peak and a suppression of the cD = ±2 peaks relative 
to the classical picture. However, in the rj = —tt/2 case they have a much more profound 
effect on the noise spectrum. Ref. |TTj shows that in this case the current noise near cD = ±1 
is proportional to the momentum spectrum of the oscillator (for our momentum-coupled 
system, this would be the position spectrum under the canonical transformation), and again 
there is a leading classical term and a quantum correction, but now the quantum correction 
can be larger than the classical term for a cold enough environment (fc^T <^ eV), producing 
a negative peak at a; = ±1. Even more interestingly, there is one more term in the noise at 
1] = — 7r/2, namely the last term in Eq. (8) of Ref. [TT], which is non-negligible when the total 
oscillator damping due to the environment and the detector is very small, as is the case in 
our plots. This is exactly the term leading to the Fano-like resonance-antiresonance features 
in our spectra. Thus, due to the non-zero tunneling phase, our spectra show clear signatures 
of the correlations between the junction current and back-action force on the oscillator. 

It is also interesting to note that our spectra are similar to those obtained classically 
by Armour for the full noise spectrum of a SET whose capacitance depends linearly on 
the position of a nearby nanomechanical oscillator.'^! In both cases, the same peaks and 
resonance-antiresonance features emerge, and the dependence of the peak heights and widths 
on the system parameters is very similar. It is possible that the quantum effects in our 
system somehow mimic the classical effects in the SET-oscillator system. The resonance- 
antiresonance features at a; = ±1 are also predicted to appear in the back-action force 
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spectrum for an oscillator coupled linearly via its position to a generic detector.'!^ In our 
case, due to rj = — vr/2, the detector current is exactly correlated with the back-action force, 
so it is not surprising that we see the same features in the current noise spectrum as well. 

Next, we focus on the regime of low external temperature, relatively low bias voltage (but 
still ^ 1, as required by the Born-Markov approximation), and comparable external and 
internal damping, i.e. fo ^ tlX^/2TT. We start with to = 0.1, V = 100, A = 0.01, f = 0.01 
and Fq = 10~^ as the central point in our parameter space, and vary each parameter around 



its central value, keeping the other parameters fixed (Figs. [7|- [lT|. In Fig. 10, the spectrum 
does not change appreciably as a function of external temperature for T < 1, hence it is 
only plotted for T = 0.01 and T > 10. In this regime, the magnitude of the a) = peak 
ranges from 9 x 10"^ to 80, the Co = ±2 peaks range from 6 x 10~^ to 40, the u = ±1 
resonances vary between 5 x 10~^ and 2 x 10~^, and the antiresonance features vary from —1 
to —2 X 10^^. The dependence of the peak heights and widths on the parameters is much 
the same as in the high-voltage, high-temperature regime. The main difference between the 
two regimes is that the cD = ±1 antiresonance features are relatively more prominent in the 
low-temperature regime. In fact, the set of parameter values to = 0.1, V = 100, A = 0.01, 
T = 0.01 and Fq = 10~^ appears to be a crossover point in parameter space, where all the 
peaks are similar in magnitude (e.g. the a) = peak magnitude is 0.3, the u = ±2 peaks are 
0.15, and the cD = ±1 antiresonances are —0.25; the resonances are negligible throughout 
this regime). For higher values of to, V, A and T, and lower values of Fq, the cD = 0, ±2 
peaks dominate, and for lower values of to, V^, A and T, and higher values of Fo, the a; = ±1 
antiresonances dominate. Besides, unlike in the high-temperature, high-voltage regime, here 
the current noise spectrum is sub-Poissonian {ASj < 0) for |a)| > 1, except near the u = ±2 
peaks, especially for large to, V or A, and approaches zero from below as cD — > ±oo. This 
is consistent with Ref. [H], where the last term in Eq. (8), which is responsible for the 
Fano-like feature at a) = ±1, grows with to, V and A. Also, the first (Lorentzian) term is 
expected to be negative when ksT <^ eV, explaining the suppression of the resonance peak 
in this regime. 

Finally, we investigate the regime of high external temperature (T = 10^) and relatively 
low bias voltage {V = 10^), and vary each of the remaining three parameters around the 
parameter-space point with coordinates to = 0.1, A = 0.01 and Fo = 10~^, keeping the 



other two parameters fixed (Figs. 12 - 14). We find good qualitative agreement between 
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FIG. 7. Dimensionless non-Poissonian current noise spectrum for to = 0.05, 0.1 and 0.2. The 
remaining parameters are V = 100, A = 0.01, T = 0.01 and Tq = 10~^. 




FIG. 8. Dimensionless non-Poissonian current noise spectrum for V = 50, 100 and 200. The 
remaining parameters are to = 0.1, A = 0.01, T = 0.01 and Tq = 10^^. 
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FIG. 9. Dimensionless non-Poissonian current noise spectrum for A = 0.005, 0.01 and 0.02. The 
remaining parameters are to = 0.1, V = 100, T = 0.01 and Tq = 10~^. 
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FIG. 10. Dimensionless non-Poissonian current noise spectrum for T = 0.01, 10 and 100. The 
remaining parameters are to = 0.1, V = 100, A = 0.01, and Fq = 10~^. 
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FIG. 11. Dimensionless non-Poissonian current noise spectrum for Fq = 10~^, 10"'', 10~^ and 
10~^. The remaining parameters are to = 0.1, V = 100, A = 0.01, and T = 0.01. 

this regime and the high-voltage, high-temperature one presented in Figs. [2]-|6} The u = 0, 
±2 peaks are again much higher and broader than the a) = ±1 resonance-antiresonance 
features, and the magnitudes of the peaks and their dependence on the parameters are very 
similar in both regimes. This agreement is not surprising given that the two regimes have 
very similar values for all parameters except the voltage. Apparently, simply lowering the 
bias voltage, while keeping it large compared to the oscillator frequency, does not introduce 
us to a fundamentally new regime. Nevertheless, there are some subtle differences. First, in 
the low-voltage regime the antiresonance features are even less pronounced and disappear 
completely for most parameter values, as seen in the figures. This is expected, as a high 
external temperature can wash out the noise suppression due to a weak (low bias voltage) 
detector, especially at weak coupling, weak tunneling and/or high external damping. Sec- 
ondly, there are slight departures from monotonicity in the dependence of the peak heights 
on the coupling and external damping parameters. The peak heights actually decrease as 
A increases from 0.01 to 0.1, and increase as Fq increases up to 10~^ before eventually 
decreasing. At present, it is not quite clear what causes these slight aberrations. 

In the experiment of Ref. [16], the parameter values (e.g. for sample A) are V ~ 10^, 
A = 7.2 X 10"^, fo = 3.3 X 10"^, and T = 1.6 x 10^,^^ corresponding to our high-voltage, 
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FIG. 12. Dimensionless non-Poissonian current noise spectrum for to = 0.001, 0.01 and 0.1. The 
remaining parameters are V = 100, A = 0.01, T = 10^ and Tq = 10~^. 




FIG. 13. Dimensionless non-Poissonian current noise spectrum for A = 0.001, 0.01 and 0.1. The 
remaining parameters are to = 0.1, V = 100, T = 10^ and Tq = 10^®. 
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FIG. 14. Dimensionless non-Poissonian current noise spectrum for f o = 10"^ 10""^ and lO'^. The 
remaining parameters are to = 0.1, V = 100, A = 0.01 and T = 10^. 

high-temperature regime, with the external damping much larger than the detector damping 
(Fo ^ tl\'^/2n). It is important to note, however, that the tunneling to is much larger 
in the experiment (tg = 0.5). Comparing the experimental and theoretical results (e.g. 
Fig. 4d in Ref. [16] and Figs. |2] - [6] in this article), we find the same peaks at a) = 0, 
±1 and ±2. Both the experimental Fano factor and the theoretical current noise vary over 
many orders of magnitude, indicating strong electron-electron correlations due to interaction 
with the oscillator. The Co = ±2 peaks are somewhat less pronounced in the experimental 
results. The small antiresonances predicted theoretically in this regime are not resolved in 
the experiment, possibly due to the background noise. The sub-Poissonian noise observed at 
higher frequencies in the experiment is absent from Figs. [2]-[6| but surprisingly does appear 
in the low-temperature, low- voltage regime in Figs. [T]- [TT| 



For comparison purposes, in Fig. 15 we have also plotted the theoretical current noise 
spectrum for the exact set of parameter values used in the experiment. The theory fails to 
predict the high super- Poissonian peaks, the resonances at a) = ±1 or the sub-Poissonian 
noise seen in the experiment. It also overestimates the relative magnitude of the u = ±2 
peaks. These discrepancies are probably due to the breakdown of the theoretical method 
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at strong tunneling, and suggest that comparisons between theory and experiment should 
be made with caution until the strong tunneling regime of our system has been investigated 
theoretically. As expected, our theory fails to predict the sharp increase in peak heights as 
the system enters the strong tunneling regime. In fact, making the coupling larger and the 
external damping smaller than in the experiment, as in Figs.|2]-|6j partially "compensates" for 
this failure, producing higher peaks more similar to those seen in the experiment. 
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m 

FIG. 15. Dimensionless non-Poissonian current noise spectrum for the parameter values used for 
sample A in the experiment of Ref. [16j: to = 0.71, V = Ix 10^, A = 7.2 x 10"^, f = 1.6 x 10^ 
and fo = 3.3 x 10"^. 

V. WIGNER FUNCTION REPRESENTATION AND STEADY-STATE OSCILLA- 
TOR DYNAMICS 

In order to investigate the steady-state behavior of the oscillator, it is useful to trans- 
late the master equation into the Wigner-Weyl formalism.^ Using the method outlined in 
Bennett,'^ one obtains: 

dtW{x, p; x; t) = I - -^d^ + mul^xdp + Fo ( - |) 

+ D+dl + 2m'ul^+d^x + Dodl + 270^^^ 
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- ^"'ft'"^^ 1^ + "''<l*(p9p + d.x)] }w(x.p; x; i). (27) 

Here we have assumed large bias voltage (the a = —1 terms are set to zero). In what follows, 
we also assume zero temperature in the leads, just as in the previous sections. Just as in 
the case of the position-coupled oscillator and QPC studied in Ref . [9] , a careful inspection 
of the above equation shows that exactly half of the back-action damping and diffusion 



(terms involving 7+ and -D+ on the second line of Eq. (27) as well as the very last terms 
on the third and fourth lines) is correlated with tunneling (multiplied by e*-^), whereas the 
other half is independent of tunneling events. The precise correlation between the electron 
tunneling events and the momentum kicks imparted to the oscillator suggests that there is a 
departure from the simple model of the detector as a thermal bath. On the other hand, the 
uncorrelated half of the back-action means that one gains information about the oscillator 
even when no electrons are tunneling. 



The first two lines of Eq. (27) represent a classical Fokker-Planck equation for an oscillator 
coupled to two equilibrium baths, and agree well with Eq. (6a) in Ref. pj, except that in our 
case the external bath is position-coupled whereas the detector bath is momentum-coupled. 
Besides, the average back-action force on the oscillator, Fo(— 7r/2), is non-zero in our case. 



The first four terms on the third line of Eq. (27) combine to give an oscillator-dependent 



tunneling rate through the QPC, analogous to the classically fluctuating tunneling rate 
represented by the last term of Eq. (6a) in Ref. [9J. However, there are subtle differences - 
instead of a quadratic dependence of the classical tunneling rate on the oscillator position, in 
our case there is a quadratic dependence on the oscillator momentum, as well as a linear term 



proportional to dxW. Finally, the first two terms on the fourth line of Eq. (27) correspond to 
the last two terms in Eq. (6b) of Ref. |9], which represent quantum corrections to the average 
tunneling rate and arise from the difference between tunneling processes involving absorption 
or emission of a phonon. There are again some differences - in the position-coupled system, 
these terms involve dxW and xd^W, whereas in our case they are proportional to xW and 
pdpW. 

To study the steady-state dynamics of the oscillator, we need to trace over N, the number 
of electrons that have tunneled through the detector, which is equivalent to setting x = in 



the above equation (cf. Eq. ([9])). The resulting simpler equation can be integrated by parts 
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to yield coupled equations for the oscillator moments, e.g. up to second order: 

-fA^) = 2m2<7+ X +-{p) = 

at h m 

— (p) = -mu;^(x) - 27o(p) = 

d . 2, hH'^eV 2htotieV , , , o o ~ , 2,, 
dt h h m 

— ixp) = -'^^^—{p) - mul^{x^) - 2{m^ul^^+ + %){xp) + —{p^) = 
dt n m 

= 2Do - 2mul{xp) - 47o(p2) = (28) 

As expected, these equations are the same as Eqs. ( [C22| ) and ( [C23| ) in Appendix |Cl In the 
above equations, {xp) represents the symmetrized moment (a;p + px)/2, which is why the 



imaginary terms from Eqs. (C22) and (C23) are absent. We will stick to this convention for 
the rest of this section. It is interesting to consider the limiting cases in which either the 
detector or the environment decouples from the oscillator. In the former case (7+, \^ — )■ 0), 
one obtains {{xp)) = 0, ((x^)) = kBT/mco'^ and ((p^)) = mkBT (when ksT ^ hi^m), 
consistent with the equipartition theorem. Note that we are working with the irreducible 
moments, i.e. the variance and covariance, hence the double brackets. In the latter case 
(70, -Do 0), one also obtains equipartition results: {{xp)) = 0, {{x^)) = kBT^et/'^n^m ^^"^ 
((P^)) = f^^BT^eti where we have defined the detector temperature Tdet = eV/2kB- 

In fact, one can solve the moment equations algebraically in the general case. One obtains 
the following results: 

2htotieV% 



{x) 
(P) 
{{x')) 
{{xp)) 



huj^{Am'^%% + 1)' 
mhtotieV 
/i(4m27o7+ + 1) ' 
2\\ _ m^hHieVji^i + 4m^a;^7o7+ + uj) + 2hDo 
Am?ujlJi{Am?^Q^+ + l)(7o + m?uj'^^+) 

m{%hHleV - 2/17+Do) 
2h{Am?%^+ + l)(7o + m?ujl^^+) ' 
2/iL'o (4771^707+ + 4m'^a;^7^ + 1) + m^ul^hHleV 



^^^^^ 4/7(4m27o7+ + l)(7o + m2a;^7+) ' ^ ^ 

The non-zero values of (x) and {p) imply that the oscillator is in a boosted frame as a result 
of our use of an effective model for the physical system of Ref . [I6] . It is interesting to know 
if the simultaneous interaction of the oscillator with the position-coupled external bath and 
the momentum-coupled detector bath can put the oscillator into a so-called squeezed state - 
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a clear signature of quantum behavior in the system. If the initial state of the oscillator has 



a Gaussian Wigner function, the quadratic form of Eq. (27) ensures that it remains Gaussian 
for all time. Such a state will obey the position-momentum Heisenberg uncertainty relation, 
i.e. in dimensionless form 

V^Vp > 1, (30) 

where 14 = ((5^)) and Vp = {{p"^)) are the variances in dimensionless units. For a squeezed 
state, we need either Vx < I or Vp < 1.1^ The values of Vx and Vp depend on the orientation 



of the X and p axes in phase space. Fig. 16 shows that Vx or Vp is minimized when the 
covariance Vxp = {{xp + px))/2 is zero. Thus we need to find a new set of axes in phase 
space such that Vxp = 0, and therefore Vx (or equivalently Vp) attains its minimum value. 
Mathematically, this is analogous to diagonalizing the variance-covariance matrix of the 
Gaussian state. If the smaller of the eigenvalues thus obtained is less than unity, we have a 
squeezed state. 

Using this criterion, we don't find any evidence of quantum squeezing in our system in 
any of the parameter regimes discussed in the previous section. The variances are larger 
than unity and approximately equal, and the covariance is much smaller than the variances, 
resulting in equal, larger-than-unity eigenvalues. The absence of quantum squeezing is not 



surprising given the form of Eq. (27). As already discussed, when we set x = 0, only the 
first two lines of the equation remain, which represent a classical Fokker-Planck equation for 
an oscillator coupled to two independent reservoirs - the position-coupled environment and 
the momentum-coupled detector. Thus an unconditional measurement of the steady state 
of the oscillator cannot be expected to yield any quantum results. Only a measurement 
conditioned on a certain history of the current can show departures from the effective bath 
model due to the correlation between the tunneling and back-action already present in the 
master equation.'^ 

A natural question to ask is whether the covariance can ever be non-negligible compared 
to the variances. To answer this question fully, we look at the analytical expressions for the 



variances and covariance, i.e. the dimensionless versions of Eq. (29): 



V 



V 



^ifX^VjATl + HfX'To + 1) + fo coth (^) 

(ftgPfo + l)(fo + ^tgP) 

focoth(i)(ftgA2fo + j,tfX^ + 1) + ^jfX'V 



(ftgA2fo + l)(fo+^tgP) 
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FIG. 16. Phase-space diagram of a squeezed state showing that one of the variances, Vx or Vp, 
is minimized when the covciricince Vxp is zero. Th.6 solid ellipses represent the Wigiier function 
contours of a squeezed state in two different orientations relative to the axes. 

Dividing these, we obtain expressions for the ratios of the covariance to the variances: 
V.r>_ l2T2f V^-coth(i) 



^ ^tlX^il + 4r§ + HlX^To) + To coth (i) 

Vp °focoth(i)(ftgA2fo + ;^t^A4 + l) + ^tgW" 

It is easy to see that for Fq < 0.1 the magnitudes of these expressions are bounded from 
above either by the external damping 2Fo or by the detector-induced damping tgA^/vr, the 
latter being a very small quantity for the physically relevant values of the coupling and bare 
tunneling parameters. The former bound is achieved when t'^X'^V /27i ^ Fq coth(l/2T). 
Thus the ratios of the covariance to the two variances can be made as large as 20% for 
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very large external damping, Tq = 0.1, low external temperature and sufficiently large bias 
voltage, tunneling amplitude and/or coupling parameter. For example, when to = 0.1, 
A = 0.1, V = 1 X 10^ T = 0.01 and fo = 0.1, one obtains 14 = 166.5, Vp = 160.1 and 
Vxp = —31.8, leading to covariance matrix eigenvalues of 131.3 and 195.3. Note that the 
bias voltage has to be very large in order to satisfy the above condition for convergence to 
the 2ro bound. In this case, the eigenvalues are large but unequal, indicating an elliptical 
Gaussian Wigner function. 

The above result is in stark contrast with the case of two position-coupled baths and 
1] = 0, where one always obtains a circular state, whose variance is determined by the 
equipartition theorem with effective temperature Tes = {'joT + 7det^det)/7eff, as discussed 
in the Introduction. It is important to note that the eigenvalues are still between the 
equipartition theorem results for the two individual baths (in this case, Vdet = 10^ for 
the detector bath and Vex_t = 0.02 for the external bath), and this is also the case for 
all other parameter values. Besides, the naive picture that the detector bath temperature 
determines the variance in momentum and the external bath temperature sets the variance 
in position is clearly wrong - in fact, in this example Vx > Vp, even though Tdet ^ T. 
The steady state of the oscillator can be characterized as a classical, thermomechanically 
squeezed state, similar to those studied in Ref. [12]. As already discussed, there are two 
sources for this thermomechanical noise squeezing - the interplay between the position and 
momentum coupled baths as well as the non-zero average back-action force due to rj ^ 0. 

Alternatively, if we allow to and A to be as large as 0.5, we can allow V to be smaller 
while keeping the above convergence condition approximately valid - e.g. when to = 0.5, 
A = 0.5, V = 10,f = 0.01 and fo = 0.1, we have 14 = 1-85, Vp = 1.81, 14p = -0.162, 
and eigenvalues equal to 1.67 and 1.99. We have made the tunnehng and coupling as large 
as possible, the bias voltage as small as possible (but still large compared to the oscillator 
frequency), the temperature as low as possible, and the damping as large as possible (but 
smaller than the oscillator frequency so that Q ^ 1). Our intuition suggests that this is as 
close to a quantum regime as we can push the parameters and still maintain the validity 
of the Born-Markov approximation. The variances in this case are close to the Heisenberg 
uncertainty limit, the covariance/variance ratios are large (about 10%) and the eigenvalues 
are unequal, yet even in this extreme regime what we have is simply a very cold, classical, 
thermomechanically squeezed state, consistent with our expectations. 
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VI. CONCLUSION 



We have studied the current noise spectrum and steady state behavior of a resonator 
coupled hnearly to a QPC via its momentum for a wide range of system parameters. Our 
spectra show clear signatures of the non-Gaussian correlations between the junction current 
and the back-action force on the oscillator, namely the resonance-antiresonance features at 
a; = ±1. These features are prominent in our case because the tunneling phase is set to a 
value (?7 = — 7r/2) where the current and back-action force are maximally correlated. Our 
results are consistent with the analysis of Ref. [11], implying that, as far as the current noise 
is concerned, the momentum-coupled system is quite similar to the position-coupled system 
with the same rj due to canonical invariance, despite the presence of a position-coupled 
external bath, at least in the case of a weakly coupled environment {Q ^ 1). 

Comparing our results to the experimental noise spectra obtained in Ref. [15], we find 
that inserting the experimental parameter values into our calculation fails to reproduce some 
important features of the experimental results, such as the high super-Poissonian values at 
the peaks or the sub-Poissonian noise away from them. This breakdown of the theory 
at strong tunneling is expected as both the polaron transformation and the Born-Markov 
approximation we have used hinge on the assumption of weak tunneling. However, if we keep 
the tunneling weak, and use stronger coupling and weaker damping than in the experiment, 
we can obtain spectra much more similar to the experimental ones, suggesting that a future 
theoretical approach that does not depend on the weak tunneling assumption might be 
much more successful in predicting the noise quantitatively. In future work, we plan to use 
scattering matrix methods to treat arbitrarily strong tunneling, as proposed by Bennett et 

aim 

Our study of the oscillator steady state indicates that once the detector is traced out, the 
oscillator obeys a classical Fokker-Planck equation, where it is coupled to two independent 
reservoirs. Thus an unconditional measurement of the steady-state oscillator moments is 
not expected to yield any deviations from classicality. However, the full master equation 



for the coupled system, Eq. (27), clearly contains quantum terms showing that exactly half 
of the electronic back-action is correlated with tunneling. Ref. [9] suggests that in order to 
observe these departures from the effective bath model, one needs to study the conditional 
evolution of the oscillator, based on a certain current measurement history, or else look at 
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the current noise spectrum of the detector, as we have done in the present study. 

Despite its classical nature, the oscillator steady state can experience significant thermo- 
mechanical noise squeezing. At high external damping, low external temperature, and large 
bias voltage, tunneling amplitude and/or coupling strength, the variances in position and 
momentum can diff'er by up to about 20%, while still remaining between the limits set by the 
temperatures of the two individual baths. The sources for this thermomechanical squeezing 
are the simultaneous presence of a position-coupled and a momentum-coupled bath, as well 
as the non-zero average back-action force on the oscillator due to rj = — 7r/2. 

To conclude, there are several future directions that one can take in order to extend the 
present study. First, one could use a scattering approach to treat the case of strong tunneling 
and enter the parameter regime of current experiments.^" Secondly, using scattering or some 
other approach, one could attempt to relax the high bias voltage {eV/hum ^ 1) assumption 
and enter a regime where one might expect to see quantum signatures in the oscillator steady 
state as well as in the current noise spectrum. Thirdly, one could look at the conditional 
evolution of the oscillator based on a certain measurement history of the QPC current, which 
can be highly non-thermal even in the weak tunneling limit .'^ 
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Appendix A: Polaronic Transformation 

Starting with Eq. ([T]), we perform a polaron-like transformation on the Hamiltonian, i.e. 
H — UHW , with the unitary operator 



e'GiOe-^Gi ^ o + O] + ^[G, [G, O]] + ■ ■ ■ + [G, . . . [G, O]] . . .], (A2) 




(Al) 



We make use of the Baker-Hausdorff lemmaP^ 
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where O is an operator, G is a Hermitian operator, and ^ is a real number. In our case, 
we identify ^ = Xxzp/{2hwm) and G = iiJ^L^^L^L — J^R^^R^Rli^'' " '^)- Expanding to first 
order in the oscillator coordinates, i.e. in ^, and using the canonical commutation (anti- 
commutation) relations for the a^^-* (b^P) operators, we obtain the following results for the 
various terms in the original Hamiltonian: 



U( Y,f^^LRb{bRW 



L,R 



1 - 



hhJr. 



(a^ - a) 



J2^^LRb[b 

R-, 

L,R 



(A3) 
(A4) 
(A5) 
(A6) 

R R 

To obtain the last two equations, we have dropped the second term in the Baker-Hausdorff 
formula as it leads to terms quartic in the b^ and bji operators. As the L and R reservoirs 



U 



U 



\x ( 

^ L R 

J2i^L - \x/2)b{bL = J2(^L - \x/2)b{bL, 

L -I L 

J2i^R + Xx/2)b^^bn] f/t = ^(e^ + Xx/2)b^^bn. 



are to be combined with the E and C reservoirs, respectively (see comment after Eq. (AlO) 
below), there is negligible accumulation of electrons in these reservoirs, hence higher order 
terms in 6^ and bn do not contribute. Also, we have 

u(Y1 ^Eb\bE + ^cb^bc) U^ = Y1 '^^EbE + Yl 'c^cbc, (AT) 

^ E C E C 

since U clearly commutes with i/bath- Finally, for the remaining two terms in -ffmt, we get 



^ E,L ' 

u(Yf^^cRblbR)u^ 



C,R 



1 H — [a' — a) 

ZTVjJm 



1 — — — —[a^ — a) 



2fUxJr. 



hriELbEbL, 

E,L 

bhbn, 



(A8) 
(A9) 



C,R 



where we can neglect the momentum-dependent second terms. 

Putting all the terms in UHW together, we see that the oscillator position-dependent 
terms cancel, and obtain 

UHU^ = tkuma^a + ^ eLb{bL + ^R^r^r + Yl ^Eb^E^E + Y ^cbcbc 

Y^^LRb\b R 



+ Y ^Eib^bL + Y ^^CRbi^bn + 

E,L C,R 

+ H.c. of last 3 terms. 



L,R 



(AlO) 
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Finally, we assume that the lead-reservoir tunneling amplitudes Qel and flcR are much 
larger than the tunneling amplitude ^Ilr between the two reservoirs, so we can effectively 
combine the emitter and the left reservoir, and also the collector and the right reservoir. 
Thus the transformed Hamiltonian can be written in the much simpler form H = iJosc + 
-f^bath + -f^int, where 

Hose = hcUmO'^O', 
L 

Ax. 



R 



Y,^^LRb{bnY + }l.c. 



(All) 



L,R 



Appendix B: Detailed Derivation of the Born-Markov Master Equation 



The starting point for the derivation of Eq. (10) is the general Born-Markov master 
equation 

■^Poscit) = T^[Ho,posc{t)] - ^ J C?t'Trbath{[^int, [^int(-i')>Posc(i) ® Pbath]]}, (Bl) 

where Hq = Hose + -f^bath and Hi^t{—t') is given in the interaction picture. The derivation 
proceeds in exactly the same way as in Eqs. (7.49)-(7.64) of Ref. [10], since the position or 
momentum dependence of the coupling is not made explicit until later. Thus we can take 
Eqs. (7.63)-(7.64) (adapted to our notation) as the starting point of our calculation: 
d 



where 



1 1 r°° 

^^Pix;t) = -^[Ho,pix;t)] - I dt'Y:Aix,R,L;t,t') 



A{x, R, L- 1, t') = [rr+(-t')p(x; t) - r^{-t')p{x; t)7le-^(^--«)*'/^;^(i _ 

+ [T^T{-t')p{x:t) - r(-t')p(x;t)rt]e^(^--«)*'/ry^(i _ 
+ [p(x; t)T\-t')T - Tpix; t)r\-t')]e-^^^^-^-^''^'fLii - Ir) 



i)[rp(x;t)rt(-t')e-^(^--^-)*'/VL(i - Jr) 



[e^^-l)[r{-t')p{x;t)'Pe 



iieL-eR)t'/h- 



- e 



i)[T^p{x;tm-t')e 

l)[r\-t')p{x;t)re 



Wl-fR) 



-i{eL-eR)t' /Tl 



]M1-/l). 



(B2) 



(B3) 
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In this equation, 6^(6/?) and fhifR) are the energy levels and Fermi finctions of the left 
(right) reservoir and T is defined as T = to + tie^^p such that 

f = ^rrt. (B4) 

Using the interaction picture expression for the harmonic oscillator momentum, 

p{t) = —mUmXsm{umt) +pcos{uJmt), (B5) 



we can write 



r(t) =to + '^-^e'\p + imuj^x)e'^-' + |e^^(p - imuj^x)e-'^-\ (B6) 

Zi Zi 



with an analogous expression for T^(t). Substituting these into Eq. (B3), we first perform 
the integration over i! using 



oo -| 







dte^''^' = Tt6{uj) ± ipv ( - . (B7) 



In the case of position-dependent coupling, the principal value term leads to a term that 
renormalizes the oscillator frequency plus another negligible term in the Caldeira-Leggett 
equation. In our case, we expect it to lead to oscillator mass renormalization (i.e. an extra 
term oc p^) plus a similar negligible term, so we drop the principal value part in complete 
analogy with the position coupling case. Next, we replace the summations over discrete 
energy levels by integrals over energy with constant density of states A: 

/■oo 

V . . . ^ / rfe^A . . . , (B8) 

where i = L,R. The result of the integration over the reservoir Fermi functions can be 
expressed in terms of the tunneling rates T±{E): 

POO 

hT^iE)= / de\to\'fie-fiL)[l-fie-fiR + E)], (B9) 
Jo 

POO 

hT4E)= / de\to\'fie-fin)[l-fie-fiL + E)], (BIO) 
where fii is the chemical potential in reservoir i. After a somewhat lengthy but straightfor- 



ward collection of terms, one arrives at Eq. ( 10 ) 
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Appendix C: Detailed Solution of the Master Equation to Find the Current Noise 
Spectrum 



In this appendix, we solve for the time dependence of the cumulants {{xN)) and {{p'^N)), 
and integrate the MacDonald formula to obtain exact analytical expressions for the non- 



Poissonian current noise spectrum ASi{uj) in Eq. (23). First we solve for {{xN)). Using the 
method of taking derivatives with respect to x then tracing the master equation over 
the oscillator degrees of freedom, we obtain two coupled first-order differential equations for 
{{xN)) and {{pN)): 

|((xiV)) = -2mW^%{{xN)) + ^{{pN)) + - mW^j^{x) 

We replace all averages that do not contain with their stationary values in order that the 
current-current correlation function {I {t + t') I (t)) be independent of t. Eliminating {{pN)) 
and imposing the boundary conditions {{xN{t = 0))) = {{pN{t = 0))) = 0, we obtain the 
following initial value problem for y = {{xN)): 

y + ay + by = c, 
y{0) = 0, 

2/(0) = rf, (C2) 

where 

a = 2(m2w^7+ + 7o), 
b = + 4771^7+70), 

2heVtoti% 2-^0 o22~~/\ 2~/\ 

_ 2jneV_q%^^^ _ 4^^^^7.7o |(^^2^ _ ^^^2) _ 2rr^|^^^^^ _ ^^^^^^^ 

h mh 

heVtoti 2 2 ~ / \ iheVt] , , 
d = m uj^l+{x) — {p) 

32 



- ?!I^|M±|((x^) - (x)^) + ^((xp^) - (xm). (C3) 

This is a simple second-order linear inhomogeneous differential equation with constant co- 
efficients, whose general solution is 

y{t) = Ae'-'' + Be''^' + j, (C4) 

b 

where 

= -^-y^b ,C5) 

are the roots of the auxiliary equation (assumed to be distinct). Applying the initial condi- 
tions, we find 

+ bd 



A 



b{ri - r2) 



B = -^^y (C6) 



We substitute this solution into Eq. (23) to obtain the first term in ASj. The resulting 
integrals converge provided that Re(ri 2) < 0, which is easily seen to be the case. The term 



dtsinujti^-j (C7) 

is not integrable, but can be evaluated by the method of Cesaro summation: 

K 

Ksmujtdt = —. (C8) 

CO 

We trust this method as it removes the discontinuity in the noise spectrum at w = and 
also makes the noise go to zero at large frequencies instead of having a zero level that's very 
large in magnitude and varies randomly with the choice of parameters. Finally, one obtains 



for the first term in Eq. (23), i.e. for the contribution from {{xN)): 

.^(lstterm)_ Sm^ul^^+to 2 \ , , ^ ^ M tnn\ 

^'^i ~ fe — e w -2"- — ^ + "2"^ — 2 ~^ t:\~i ■ y^^) 

h ti lr( + u^ r2+UJ^ b\ujJ\ 
In a similar way we can solve for {{p'^N)), needed to calculate the second term in the 
non-Poissonian current noise. We first obtain three coupled linear differential equations for 
the cumulants that are second-order in the oscillator variables: 

^({x'N)) = c,{{x'N)) + C2{{{xp)N)) + Cs{{xN)) + c„ 
dt 

U{{xp)N)) = c,{{x^N)) + Ce{{{xp)N)) + c,{{p^N)) + Cs{{pN)) + cg, 
dt 

j^iip'N)) = c,o{{{xp)N)) + Cuiip'N)) + C12, (CIO) 
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where we have defined the constants 



2 
m 

2htotieV 

2htotieV o 2 2 ~ / 2\ 2ihtleV 
C4 = (x) - 2m w„7+(a; ) {xp) 



C5 



h U h 

.2 



C6 = -2(m2a;^7+ + 7o), 
1 

C7 = — , 

m 

htotieV 

cg = ihm'^u'^^+ + im^ul^%Y^{x) + ^^^^ {p) - 2m^ul^^+{xp) 

ihtjeV , n, 2m^w£,7+to/, o \ / \ / w ^i^^r, s\ / \/ 2\\ 



Cio = —2muj'^ 



m' 



h h ti h 

rr 

cii = -47o, 

ci2 = 2.m^c.^7+^(p) - 2mW^%{p') - ^^^^^^((xp^) - (x)(/)) 
ti n t\ 

+ ^((/) - (Cll) 



Ehminating {{x'^N)) and {{{xp)N)) and substituting for {{xN)) and {{pN)) using Eqs. (C4) 
and ( |C1[ ), we get a third-order hnear inhomogeneous ordinary difi^erential equation for z = 
{{p'N}): 

where we have defined the constants 

a = -(ci + C6 + Cll), 
/3 = CiCg + CiCii — C2C5 + CgCii — C7C10, 
7 = -cicecii + C1C7C10 + C2C5C11, 
H = ^(cgCsCio + CsCiomrl + 208010171^^1^^ +ri 
- CiCgCiomri - 2ciC8Ciom^u;^7+), 



34 



V = S(c3C5Cio + c^cwmr^ + 2c^ciQmru^^+r2 

- CiCgCiomra - 2ciC8Ciom^a;^7+), 

p = c^c^ciQ- - 2ciC8Ciom^a;^7+- + ciC6Ci2 + ciCgCiomd 

- C1C9C10 - C2C5C12 + C4C5C10, (C13) 
and b, c, ri, r2, ^ and B are as defined above. The general solution is easily found to be 

z{t) = Ce^'i* + DeP^^ + Ee"^^ + Me'''* + A^e''^* + ^ , (C14) 

7 

where Pi are the three roots of the auxiliary equation + ap^ + /3p + 7 = (assumed to be 
distinct), the constants C, D and £^ are to be found by imposing initial conditions, and 



r\ + arf + /3ri + 7 ' 



1/ 



2 5 • (C15) 

+ ar^ + /3r2 + 7 



Assuming — 3/3 7^ and using the cubic equation formula, we find 

a R -?>(5 



Pi 



3 3 3i? 



a 



^'3^6^ 6/? 

a , i?(l-^v^) , (l + iV3)(a^-3^) 
^^ = -3+ 6 + 6^ ' ^^'^^ 



where 



Q = ^y{2a^ - 9a/3 + 277)2 - 4(^2 _ 3/5)3^ 

= y^(g + 2Q;3- 90:^ + 277), (C17) 

and any determination of the complex-valued square and cube roots can be used. From the 
boundary conditions ((a;2A^)) = {{{xp)N)) = ((p^iV)) = {{xN)) = {{pN)) = at t = we 
obtain the appropriate initial conditions: 

z{0) = 0, 
i(0) = S, 

z(0) = e, (C18) 

where 

S = C12, 

e = C9C10 + C11C12. (C19) 
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Applying these initial conditions and solving the resulting system of equations, we obtain 
C 

D 

E 

{P2 - P3){P3 - Pi) 

(C20) 



A/7 ( /l,^ /If, 

ivi \p2P3 


' 1P2 


'iPs + ' l) 


-t- jv \p2P3 


' 2P2 


' 2P3 + ' 2) 


^ PP2P3II ^ Op2 -f- Ops t 








(Pi- 


P2)(P3 


-Pi) 






- ripi 


- rip3 + rf) 


+ A^(piP3 


- ?^2Pl ■ 


- r2pz + rf) 


+ PP1P3/7 + ^Pi + Sp3-e 








(Pl- 


P2)(P2 


-P3) 






- ripi 


- rip2 + rj) 


+ A^(PlP2 


- '^2Pl - 


- r2P2 + A) 


+ PP1P2/7 + ^Pi + 5p2-e 








(P2 - 


P3)(P3 


-Pi) 





Having determined all the relevant constants, we finally substitute Eq. (C14) into the Mac 



Donald formula, Eq. (23), to obtain the second term in AS"/. The integration proceeds in 
exactly the same way as in the calculation of the first term. One has to check that the roots 
Pi, P2 and ps have negative real parts, which is best done numerically as the analytical ex- 



pressions are rather complicated. Finally, we obtain for the second term in Eq. (23), i.e. the 
contribution from {{p'^N)): 

. ^(2nd term) _ AeW^ [ Cu Du Eu Mu Nu p / 1 

- ^ ^2 , , ,2 + ^2 I , ,2 + ^2 I , ,2 + ^2 I , ,2 + ^2 I , ,2 + ~ I ~ 



(C21) 



In the above calculations, we needed to know the A^-independent oscillator moments, 
), up to fourth order. These can be found easily by solving systems of linear equations 
successively, starting with the first-order moments, then going to second order, etc. The 



equations for the moments are obtained by miltiplying Eq. (10) by x^p' , setting x = 0, 
tracing over the oscillator degrees of freedom, and finally setting d{x^p')/dt = since all 
oscillator moments are stationary as already discussed. The moment equations are given 
below. 

First-order moments: 

d , , HtoticV 2 2 ~ I \ / \ 

37W = r 2m w^7+ a; + — P = 

dt h m 

j^{p) = -mul{x) - 2%{p) = (C22) 



Second-order moments: 

c? . 2\ hH'^eV 2htotieV .22~/2\ 2,, 

-ix^) = + + \ {x) - Am^ul^+{x^) + -{xp) = 

dt m h n m 

— {xp) = ih{m'^ul^^+ + 7o) + (p) - mul^ix^) - 2(m^a;^7+ + %){xp) + — (p^) = 

dt n m 

j^) = 2D, + thmul - 2mujl{xp) - A^l^p") = (C23) 
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Third-order moments: 

d , o, /lT?t]eV ih\ , , ^htntieV ,9, 99 ,0, 3 , , 
-{x') = 3(— i + y (x^) - 6m2^^7+(^') + -(^^'p) = 

— {xp}= 2in[m a;^7+ + 70) (a;) + — — j (p) + {xp) 

777- 

^(V) = (^/imcu^ + 2Do)(a;) + 2ih{m^ujlfi+ + 70) (p) + ^^^^^(j,^) 
- 2mujl{x''p) - 2{m^ula^ + 27o)(V) + -(p') = 

= Ziihrnul + 2Do)(p) - "imwUxp") - 67o(p') = (C24) 



Fourth-order moments: 

at \ a m/ n 

- miulix^) - 2{3mW^% + ^o){x''p) + -(xV) = 

j^{xV) = {ihrnul + 2Do){x^) + m{mW^% + %){xp) + (^^^ - ^) (p') 
+ ^^^^^(V) - 2mujl^{x^p) - A{m'^ujl^% + 70) (^V) + —{xp^) = 

ft TDj 

j^{xp^) = 3{ihmujl + 2Do){xp) + 3ih{mW^% + %){p^) + ^^"^^^^^ (^3^ 

- Smulixy) - 2{m^uji% + 3%){xp^) + -{p'') = 

^(p^) = 6{thmujl + 2Do){p^) - 4ma;^(xp'^) - 87o(p'^) = (C25) 



Importantly, the equations for the moments of a given order only involve moments of lower 
order, which have already been calculated and become part of the constant vector in the 
resulting matrix equation. Thus there is no need for truncation or use of a semi-classical 
approximation in our case. 
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